The Role of Baryons in Unified Dark Matter Models 
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We discuss the importance of including baryons in analyses of unified dark matter scenarios, fo- 
cusing on toy models involving a generalized Chaplygin gas. We determine observational constraints 
on this unified dark matter scenario coming from large scale structure, type la Supernovae and CMB 
data showing how this component can bring about a different behaviour from classical ACDM and 
thus motivate further studies of this type of models. We also speculate on interesting new features 
which are likely to be important on non-linear scales in this context. 



I. INTRODUCTION 

Unified dark matter (UDM) models — also referred to 
as 'quartessence' — where dark matter and dark energy 
are seen as different manifestations of a single substance 
have recently enjoyed considerable attention, the main 
reason for it being, to try and find viable alternatives 
to the highly fine-tuned potentials so characteristic of 
quintessence models. The price one has to pay seems to 
be the introduction of fluids obeying exotic equations of 
state. One such example is the class of models based 
on the Chaplygin gas, first suggested for this purpose in 
PI Q and trivially generalized in Q . 

There has been some controversy about the viability of 
these scenarios. Following an earlier, simpler analysis 
we were able to obtain constraints on this type of model 
from a combined dataset of 92 high- z type la supernovae 
(including forecasts for the proposed SNAP experiment) 
as well as the matter power spectrum [Bj]. This has been 
subsequently repeated by a number of authors 0,0 , both 
of which confirm our earlier results, while @ have done a 
Bayesian analysis of a sub-sample of the supernovae data 
set. The main message of these studies is the fact that 
constraints critically depend on whether one treats the 
Chaplygin gas as true quartessence (replacing both dark 
matter and dark energy) or if one allows it to coexist with 
a 'normal' dark matter component (which could be called 
the 'Chaplygin quintessence' scenario). As one might 
have expected ah initio, the case for the Chaplygin gas 
is stronger in the former scenario. 

However, the above analyses essentially restrict them- 
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selves to the background cosmology. Going beyond this 
and studying perturbation theory in this context seem- 
ingly dealt the first blow to the UDM possibility. The 
authors of [j| have studied the implications of a sin- 
gle generalized Chaplygin component to the mass power 
spectrum and by comparing it to 2dF observations were 
able to constrain parameters controlling it to such values 
as to render it virtually indistinguishable from a usual 
ACDM scenario. The production of violent oscillations 
causing the exponential blowup of the power spectrum, 
totally inconsistent with observations was the reported 
fatal flaw. This was further claimed to be of a general 
nature within UDM scenarios, which if true could be con- 
strued as evidence for independent origins of dark matter 
and dark energy. 

The main goal of this Letter is to show that such flaw 
is in fact due to ;9j overlooking a most important ingre- 
dient describing the average Universe: baryons. (This 
concern has also been voiced, but not exploited, by 0). 
Our main objection can be brought down to the follow- 
ing point: The main property (cosmologically speaking) 
of a Chaplygin gas is that of mimicking cold dark mat- 
ter (CDM) at early times as it progressively evolves into 
a cosmological constant. This results in perturbations 
in the Chaplygin gas component being heavily damped 
at late times. However, if to it we add an independent 
component with a low sound speed velocity (as baryons 
have), the normal growth of inhomogeneities can still go 
on when the Chaplygin gas starts behaving differently 
from CDM. By just considering the Chaplygin fluid the 
authors of [9j have artificially constrained the possible 
power spectra that UDM scenarios can cover. Although 
baryons are not that important for other studies such as 
constraining these scenarios from supernova luminosity 
distances, they are of critical importance in the context 
of large scale structure studies. 

We should also mention two recent papers which 
accurately study perturbation growth in these models 
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[Tol Hll | (including the cosmic microwave background). 
Although of a broader scope than ours (they consider 
a frarj/on+CDM+Chaplygin universe), we do not quite 
agree on the latter's interpretation of UDM. While 
we agree with their conclusion that this 'Chaplygin 
quintessence' scenario is all but ruled out (or at least dis- 
favored by the current data relative to ACDM) it seems 
to miss the point that these scenarios came into existence 
as an attempt to unify dark energy and dark matter, and 
that in this 'quartessence' context its behavior can be dif- 
ferent from ACDM. Yet this quartessence scenario is all 
but ignored in their discussion. Therefore, we re-derived 
the analysis of (in order to accommodate for baryons) 
and part of the analysis of ^l| to explore more fully the 
viability of the 'quartessence' scenario. 



II. GROWTH OF PERTURBATIONS 

To study the role of baryons in a UDM model based 
on the Chaplygin gas we obviously need relativistic equa- 
tions governing the evolution of perturbations in a two 
fluid system. For a general case of n gravitationally in- 
teracting fluids, the linear evolution of perturbations in 
the comoving synchronous gauge is given by |12j : 

h" + (2 + £)ti + 3 V (1 + 3vf)nA = (1) 
*■ + (! + u>i)(6i/dR + h'/2) + 3(u? - Ui)Si = (2) 



% + (1 - 3^)0, + 



<zH(l 



V 2 ^ = (3) 



where ' = d/dx, x = In a (note that |12| uses confor- 
mal time instead), h is the trace of the perturbation to 
the Friedmann-Robertson- Walker (FRW) metric, H is the 
Hubble parameter, £ = H'/H, 6i is the density contrast of 
the ith-fluid obeying pi — uiipi with an adiabatic sound 
speed Vi and a Qi element velocity divergence. Note that 
and J3J apply for all i = 1, . . . , n. 
The first noticeable point is that for baryons (which we 



are treating as an ordinary CDM fluid with aj b 



0) 



Eq. immediately renders Of, — 0bo/a. Therefore, if 
9bo = then 6 b = at all times. Putting this into Eq. 
we find that h! — —2S' b . Taking this into consideration, 
the Fourier modes of a two fluid model of baryons and a 
Chaplygin fluid with equation of state p = —Cj p a (here 
C is a positive constant and a > 0) turn out to be: 



S' b ' + (2 + £)6' b - 3/2 [fl b S b + (1 - 3aio cg )fl cg S cg } = (4) 
5' cg + (1 + Lu cg ) [9 cg /aK - 5' b ] - 3^ cff (l + a)S cg = (5) 

6> cg + (l + 3aw cg )6 cg + aUcak2 



aH(l + Lu cg ) 



6 cg = (6) 



and 



where we have used the fact that vf, g — ~auj cg 
V = — fc 2 . Thus we have three equations for three un- 
knowns: Sb,5 cg ,9 cg . Given uj cg and H (and £, fib, Qcg) 
as functions of x we can easily transform this set into 



four first order differential equations and integrate it us- 
ing a standard Runge-Kutta method. Since in the linear 
regime and deep into the matter era 6b, C g oc a implying 
S' b cg oc a, normalized initial conditions [S b , S b , S cg , cg ]o = 
[1,1,1,0] were used. As in we evolve our system for a 
plane Universe from z = 100 until today and obtain cor- 
responding transfer function But unlike in the case 
of |9(, our transfer function now comes from the bary- 
onic component and not from the Chaplygin component, 
thus avoiding in the baryonic component the violent os- 
cillations in the Chaplygin gas that were so sensitive to 
a. 

Note that a (CDM processed) scale invariant Harrison- 
Zel'dovich spectrum emerging from equality is accurately 
given by .13] 



\h\=Ak 




-1/2 



(7) 



Q = k/Th, where T — Q m oh is the shape param- 
eter (f2^j being the CDM energy fraction today), A, 
a normalization constant, [fc] = Mpc -1 and e = 
[2.34, 3.89, 16.1, 5.46, 6.71]. As previously discussed in [| 
the presence of a Chaplygin fluid (still firmly behaving 
as CDM at this epoch) only affects the shape parameter 
of the power spectrum according to T = S]^/! where 



(8) 



a and A being the parameters controlling the (general- 
ized) Chaplygin gas. Unfortunately, this simple result 
relies on the assumption of the absence of baryons. In 
order to take them explicitly into account, Sugiyama's 
shape correction 0] must be used instead 



f> h0 (l 



(9) 



where flbo stands for the baryon energy content today. 
(Note that in the absence of baryons r* = T as it should). 
Thus, we have taken Q as our initial power spectrum 
with a r* shape parameter, further processing it through 
Tfc so as to compare it to 2dF results. 



III. THE ANALYSIS 

By performing a likelihood analysis of our model of 
baryons and a Chaplygin fluid using the 2dF mass power 
spectrum we were able to constrain the (a, A) param- 
eter space assuming only reasonable priors coming from 
WMAP [13, namely Cl° b = 0.044 and h = 0.71. Note that 
at the time of recombination, the Chaplygin gas would 
still firmly behave as CDM. Therefore, standard small 
scale CMB results are to be expected when one identi- 
fies f2 m with f^. The results could conceivably differ on 
very large scales, though here they would be competing 
against cosmic variance. 
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FIG. 1: 68%, 95% and 99% likelihood contours in the (a, A) 
parameter space for a model of baryons plus a (generalized) 
Chaplygin gas, coming from the 2dF mass power spectrum. 
Note the minute ACDM region near a ~ (see text). The 
zone inside the solid lines corresponds to V* = 0.2 ±0.03 fl6]| . 
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FIG. 3: 68%, 95% and 99% likelihood contours resulting from 
a joint analysis of large-scale structure (2dF) and type la su- 
pernovae for a model of baryons plus a (generalized) Chap- 
lygin gas. The zone inside the solid lines corresponds to 
F* = 0.2 ±0.03 E3. 
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FIG. 2: 68%, 95% and 99% likelihood contours for a model of 
baryons plus a (generalized) Chaplygin gas, coming from the 
92 high-z type la supernovae. The zone inside the solid lines 
corresponds to F* = 0.2 ± 0.03 0. 

As in any 2dF data over k > 0.3 hMpc^ 1 was dis- 
carded so as to stay firmly grounded in the linear regime 
where our analysis holds. Specifically, we have evaluated 
a 500 x 10_0 x 100 data grid for A, A and a with < a < 1 
and < A < 1, where corresponding probabilities were 
found and posteriorly summed over A. Resulting con- 
fidence levels are shown in Fig. ^ where two disjoint 
regions can be found: one prominent, the other small 
(around a ~ and A ~ 0.75). We have also displayed 
the region of parameter space corresponding to a value 
of the shape parameter of T* = 0.2 ± 0.03, as per [l6| . 

The significance of this smaller area is as follows. One 
of the most noted proprieties of the Chaplygin gas class 



of models is that of reproducing to all orders a classi- 
cal ACDM scenario when a = (with 1 — A being the 
equivalent CDM fraction r2* n0 ), which we know to be in 
good agreement with observations, though still lacking a 
sound fundamental physical motivation. In our case of 
baryons plus a Chaplygin fluid, we are basically in the 
same ACDM scenario, just with a slightly higher mat- 
ter content because of them. Therefore, this small area 
around a ~ is just the ACDM component of the model 
manifesting itself. On the other hand, we have an en- 
tirely disjoint region at very high confidence level. So, 
in fact, we have provided explicit evidence of a (general- 
ized) Chaplygin gas not having to behave as ACDM in 
order to reproduce 2dF large-scale structure data. 

Now, let us re-visit high-z type la supernovae. As pre- 
viously emphasized, these are only sensitive to the model 
background evolution and not to the clustering details. 
So little information is to be gained for a from a like- 
lihood analysis of our model of baryons plus a Chap- 
lygin fluid for the supernova data on their own. In 
fact, this analysis restricts mainly the value of A since 
it more directly relates to the energy fuelingthe back- 
ground expansion — a point already made in p|. Fig. [21 
shows the result of a likelihood analysis using a combined 
dataset of 92 high-z supernova observations 0, 0, 
(see Hill 

for more details). 

However, notice that the degeneracy contours (large- 
scale structure and supernovae) of the two observables, 
are not parallel. This is of course to be expected since 
they are rather different in nature and effectively arise 
from very different redshifts. Therefore a joint analysis of 
Supernovae constraints with those of large-scale structure 
can moderately improve the constraints on a, as shown 
in Fig. |3 
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IV. CONCLUSION 

We have shown that baryons play a crucial role in the 
context of unified dark matter models involving a (gen- 
eralized) Chaplygin gas, their presence enabling them to 
reproduce the 2dF mass power spectrum without having 
to behave as classical ACDM scenarios The reason is that 
baryons can carry over gravitational clustering when the 
Chaplygin fluid starts behaving differently from CDM. 
This is of great importance since it allows the (gener- 
alized) Chaplygin gas to be a conceivable quartessence 
candidate. 

On the other hand, it is easy to see why the Chaplygin 
quintessence scenario will be observationally uninterest- 
ing, quite apart from theoretical considerations: since 
in the standard scenario there is currently no preference 
of the data for quintessence over a plain cosmological 
constant, it is clear that if we add a Chaplygin fluid to 
baryons and CDM then the available data will prefer it 
to behave as a cosmological constant as well — which is 
the result of 

Hence the present results, combined with those of 
strengthen our previous conclusions [f| . If by an indepen- 
dent method we determine the total matter density of the 
universe to be O m ~ 0.3 then in the context of this model 
we would in fact require a cosmological constant so as to 
account for the current observational results. Conversely, 
the quartessence case, where baryons are the only mat- 
ter component present, is such that its differences with 
respect to the standard case are maximal, and indeed 
in this case the A-limit is already strongly disfavored by 



observational data. 

We finally wish to emphasize that this study, as well 
as the rest of the published work so far, has been re- 
stricted to the linear regime. Interesting new features of 
this scenario may also appear on non-linear scales. An 
example would be the possibility of avoiding the cuspy 
dark matter halo profiles expected in the context of CDM 
models. We also anticipate that galaxy and cluster evo- 
lution could be modified in the context of UDM models 
based on the Chaplygin gas. This again results from the 
strong dependence of the properties of the Chaplygin gas 
on the background density. 

Moreover, in the non-linear regime the density is much 
higher than the background density, and therefore one 
can not simply infer the behavior of the Chaplygin fluid 
through the usual sound speed (which uses the back- 
ground density and pressure). So in a collapsed region 
the Chaplygin gas could still behave as CDM even though 
its background behavior may already be different. On the 
other hand, since in this case the pertinent scales are very 
small, even a tiny sound speed may be sufficient to halt 
collapse. Clearly these issues must be addressed before 
one can reach a final verdict on this scenario. We shall 
return to them in subsequent work. 
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